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Abstract. We have studied the appearance of chaos in the many-body spectrum of inter- 
acting Fermions. The coupling of a single state to the Fermi sea is considered. This state 
is coupled to a hierarchy of states corresponding to one or several particle-hole excitations. 
We have considered various couplings between two successive generations of this hierarchy 
and determined under which conditions this coupling can lead to Wigner-Dyson correlations. 
We have found that the cross-over from a Poisson to a Wigner distribution is characterized 
not only by the ratio V/Ac, but also by the ratio Ac/5. V is the typical interaction matrix 
element, 5 is the energy distance between many-body states and Ac is the distance between 
many-body states coupled by the interaction. 
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I. Introduction 

Random Matrix Theory has been successful in describing the many-body spectrum 
of interacting electronic systems An issue addressed recently is to understand 

how the transition from a Poisson to a Wigner-Dyson (WD) statistics occurs when 
interaction is switched on and what is the driving parameter for this transition . 
The purpose of this work is to describe how the interaction drives the appearance 
of chaos, i.e. how the many-body energy levels can present WD correlations. We 
shall mainly study the distribution P{s) of nearest spacings s between many-body 
levels. For WD correlations, it is well described by the Wigner surmise: P(s) = 
(7r/2s)e~'^/'*^ , where s is the level spacing normalized to its mean value. 
As a frame of reference, we will consider the recent work of Altshuler et al which 
considers the structure of the many-body states when a single particle state interacts 
with a Fermi sea. Due to the interaction, this single particle state decays by the 
creation of an electron-hole pair. The final state will be called a 3-particle state 
(2 electrons and 1 hole). Because of the interaction, the initial state has a finite 
lifetime r. In an infinite clean Fermi liquid, it is well known that the half-width of 
the state T = ?7,/2t is proportional to /Ep where e is the distance to the Fermi 
energy Ep. In a diffusive system, the effective interaction is increased since two 
diffusing quasiparticles have an enhanced probability to interact several times . As 
a result, the quasiparticle width is proportional to EcY^'^ where A is the mean 
level spacing between quasiparticle energies and E^ is the Thouless energy. 
In a finite size system, an interesting new behavior arises when the quasiparticle width 
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becomes smaller than the level spacing A. This happens when e < Ec- In this case, 
individual particle peaks can be resolved and have been observed experimentally 
[ pO[ . In this regime of non overlapping resonances, the width becomes dimension 
independent and varies as A(e/i?c)^ [§■ This description stays valid as long as the 
width is larger than the spacing between final states A3 = 4A^/e^, that is e > 
e* cx ^/EcA, so that the Fermi golden rule is applicable. In this regime, the final 
states, consisting of 3-particle states, are themselves unstable and can decay into a 
hierarchy of 5, ... , 2n+l quasiparticle states, also called states of generation {2n + l). 
This hierarchical structure of the Fock space has been discussed by Altshuler et al. 
who found interesting properties of delocalization in this tree-like Fock space. One 
important issue concerns the ergodicity of the many-particle states and their spectral 
statistics which is claimed not to be described by the WD approach (We will not 
consider here the regime e < e* where the decay is not described by Fermi golden 
rule) . This work has also been reconsidered in more recent papers ||ll|;|l2| . 




Motivated by this work, we have studied the statistical properties of this many-body 
spectrum constituted of hierarchical states coupled by the interaction matrix elements. 
We assume that in the absence of coupling, the levels obey Poisson statistics. This 
would be the case for the single particle levels (1-particle states) in a clean non 
chaotic cavity. Even in the presence of disorder or in a ballistic chaotic cavity, the 
3-particle states are described by a Poisson distribution: indeed, 3-particle states are 
formed by addition of uncorrelated 1-particle states, since they have quite different 
eigenenergies: 3-particle states and more generally (2n-|- l)-particle states thus follow 
Poisson statistics. The main goal of this paper is to describe how level correlations 
evolve from Poisson to WD as the interaction increases. 

In order to describe how these correlations set in from the coupling between states of 
the different generations, we have mainly studied the coupling between two successive 
generations. First, we have considered the coupling between one state (of generation 
1) and a dense ensemble of energy levels (generation 3). This is the well-known Bohr- 
Mottelson problem |Q . When the interaction is switched on, the spacing distribution 
P{s) deviates from a Poisson statistics and in the limit of infinite coupling, it tends to 
a limiting distribution which is intermediate between Poisson and Wigner statistics. 
We study numerically how this distribution depends on the type of coupling. This 
is done in section II. In section III, we consider the case where several states of 
generation 1 are coupled to states of the generation 3, considering that the hierarchy 
stops at this generation, which is physically relevant when the energy resolution is 
limited at generation 3 by some inelastic broadening. In this case, there is a cross- 
over from Poisson to Wigner statistics which is driven by the number of intruders 
as well as the strength of the interaction. For a large number of intruders, the WD 
correlations appear when the typical interaction matrix element becomes of the order 
of the spacing between final states. In section IV, we describe the spectral function 
(LDOS) of an intruder state. 

The case of generation 3 coupled to generation 5, or more generally the case of gen- 
eration (2n — 1) coupled to generation (2n -I- 1) (with n > 2) is more subtle. It is 
considered in section V. In this case, due to the two-body nature of the interaction, 
the states of the generation (2n — 1) are only coupled to a small number of states of 
generation (2n -1-1), so that the number of final states connected by the interaction 
is much smaller than the total number of final states. It has been argued that in this 
case the cross-over to a WD statistics should be uniquely dependent of the ratio C//Ac 
where Ac is the distance between final states connected by the interaction [^^. We 
show that it is not true: the cross-over also depends on the density of final states. 
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II. Coupling between one level and a background 

We first study the statistics of energy levels when one state (also called the intruder 
or first generation state 0) is coupled to a dense ensemble of energy levels. This 
ensemble will be also called the background or quasi-continuum, or generation 3 since 
in ref. 0, it consists of 3 quasiparticles states. Their level spacing is written A3. The 
typical coupling strength is written U . We will be interested in the regime validity of 
Fermi golden rule where U > A3 or equivalently F = ttU'^ / A3 > A3. Since we want 
to address this problem on very general grounds, we will not refer more precisely to 
the problem studied in ref. . We assume that the states within a given generation 
are not coupled directly (see comment in section V). It is known that when the 
background is described b y a WD statistics, a Gaussian coupling to an intruder does 



not change this statistics 1 13 15[| . The levels are shifted and the correlation between 
old and new levels has recently been studied by Aleiner and Matveev Here, we 
want to see how the interaction with the extra level can induce correlations between 
levels of the background, starting with an original Poissonian sequence. 
Consider one state |A) , coupled to the background {|fc)} of states obeying Poisson 
statistics via the Hamiltonian H: 

Ho = exc\cx + ^ efcc[,Cfc + ^ {Vxkc{ck + h.c.) (1) 

fc=l,JV k=l,N 

The "local" Green's function G^^{E) = {\{\E -H + iOy^\\) is given by: 

The imaginary part of this function defines the Local Density of States (LDOS) 

px{E)^Y.\(X\n)\H{E-E,,) = --lmG^,{E) (3) 

l»> 

where \n) are the exact eigenstates, with energies i?„. This LDOS, also called the 
strength function, describes the projection of the initial state |A) on the eigenstates 
\n). When the background is indeed a continuum or when the energy levels are 
broadened by some mechanism, the imaginary part Tr'^j^V\kVk\S{E — e^) can be 
replaced by: F — ttU^/A^^ where U = \/Jy^ is the typical matrix element of the 
interaction and = I/A3 is the DOS of the background. Then the LDOS has the 
Lorentzian shape: 

However, when the states of the background are discrete, this Lorentzian is just 
the envelope of a finer structure. This fine structure is obtained by looking for the 
eigenstates \n) of the perturbed Hamiltonian, which are the solutions of: 

^--SS^ (5) 

fc 

The right hand side of the equation, which is a function of E, diverges at each e^: in 
each interval [ek,ek+i], there is exactly one eigenstate £"„ (there are two more levels 
with energies E < ei and E > cn which tend to ±00 with U). 
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The levels near the intruder are more perturbed than those far from it. As a 
result the new statistics depends on the position in the spectrum. To avoid this 
complication, we first study the spectrum in the case of infinite coupling U 3> A3 so 
that the spectrum is uniformly perturbed. The eigenstates are then the solutions of: 



E 



E~ek 







(6) 



We first consider the case of a constant coupling. The distribution P(s) is shown on fig. 
1. It is intermediate between the Poisson and Wigner statistics. The level repulsion 
occurs because each of the new levels is locked between two levels of the original 
sequence {ck}- This distribution has been recently studied analytically by Bogomolny 
et al. jl^ who calculated the slope at small separation s: P{s) ~ (7r\/3/2)s. 
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Fig. 1. P(s) for a Poisson background 
perturbed by a constant, Gaussian or flat 
coupling to one intruder. 



Fig. 2. 77 vs the distance EjV to the 
intruder, in the case of Gaussian coupling. 
The thin line is the LDOS of the intruder 
(in arb. u.) over the same energy interval. 

To characterize this and other distributions throughout the paper, we calculate the ra- 
tio 7? = (/;" P(s)ds-/;° Pwd{s))I{11° Ppis)ds-J^' Pwd{s)) where so = 0.4729 is 
the first intersection point of the Wigner {Pwd{s)) and Poisson (Pp(s)) distributions 
77 interpolates between 1 (Poisson) and (Wigner), see table I. For the sake of 
comparing data, we also introduce another distribution obtained if the sum (^ would 
contain only pairs of neighboring levels and ek+i- In this case, the new eigenvalues 
would simply be given by (et + €k+i)/2. Since the {e^} follow a Poisson distribution, 
the new distribution is: P{s) = 4se~^^. This distribution, called semi-Poisson, is also 
the distribution of nearest spacings for a plasma model with short range logarithmic 
interaction 
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We have also studied the level distribution when the interaction V\k is not constant. 
We took two other probability distributions: i) a uniform (flat) distribution over a 
finite interval, ii) a Gaussian distribution. We find that the repulsion is stronger with 
a constant coupling than for flat or Gaussian couplin gs. A similar effect occurs in 
the case of an intruder coupled to a GOE background [|l3[ , in which case a Gaussian 
coupling does not affect the level statistics significantly while a constant coupling 
induces a quartic repulsion. The reason for this is the following: the eigenenergies 
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given by eq.(^) are trapped between the original levels {cfc}. With a constant coupling, 
a trapped level is repelled by two original levels with roughly the same strength. On 
the contrary, with a random coupling, two eigenenergies are closer to each other when 
the coupling is small and the repulsion is weaker. 

When the coupling is finite, the spreading width T is finite and the statistics of the 
eigenstates depends on the energy distance E to the intruder. At £' 3> F the levels are 
weakly perturbed and their statistics remains close to Poissonian. On the opposite, 
when E stays smaller than F the statistics is the same as the one obtained above 
when U is infinite, fig.(||). 

The nearest spacing distribution describes only short range level correlations. To get 
some information on longer range correlations, we have also calculated the number 
variance: S^(£') = {N'^{E)) — {N{E))'^, which measures the fluctuation of the number 
of levels N{E) in a band of width E. Contrary to the WD case which is characterized 
by small fluctuations cx \n(E), here shows large fluctuations close to 

those of a Poisson distribution (fig. ^, solid thick line): for both constant and Gaussian 
couplings, S^(£^) behaves like E—1 at large E. Since the new states are locked between 
the initial states of the Poisson sequence, their compressibility is the same as for the 
Poisson sequence, so that x = lim^j^oo ^^{E)/E — + 1. 

To summarize this section, the coupling of a background with one extra level, even 
when infinite, induces a cross-over from a Poisson to an intermediate statistics which 
is still far from Wigner distribution. This is because the sequence of new levels 
alternates with levels of a Poisson sequence. 



III. Coupling between several levels and a background 



We now study the case where several levels are coupled to the background. Consider 
a set of m intruders {Ai, A2, . . . , A™} (m ^ N; N is the number of states in the 
background) with mean level spacing A. As in the case of one intruder, the Green's 
function can be calculated exactly. It is the solution of the linear m x m system: 



[{EI -Hx-M] 



( %>^^ \ 

Ai A2 



\ 


f M 







) 





with Mij — 



Vx-kVkXj 
E~ek + iO 



(7) 



H\ is a diagonal matrix, whose diagonal elements are e^i , • ■ • , • I is the identity 
matrix. The new eigenstates are given by the equation: det[{EI — H\ — M] = 0. The 
problem has thus been reduced to the diagonalization of a m x m matrix. 
Each intruder is broadened into a Lorentzian. It is clear that, as long as the width 
F of each Lorentzian is small compared to the distance A between intruders, the 
overlap between them can be neglected and we are back to the previous problem of 
a single intruder. This is the case as long as F < A, i.e. when the typical interaction 
U < y/AAs (In the case of the disordered Fermi liquid where U ~ A'^/Ec §,@,|l8| and 
A3 ~ AVe^ this means e < E^ M). 
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Fig. 3. r] vs E/A, when 5 intruders 
are coupled to a Poisson background by 
a Gaussian coupling. The energy inter- 
val is restricted between the 2"'' and the 
3'^'* intruder. Inset: Total spectral func- 



5 10 15 20 

number of intruders 

Fig. 4. Variation of rj with the number 
of intruders in the limit of infinite (Gaus- 
sian or uniform) coupling. The inset shows 
P{s) for the case of 20 intruders, compared 
with the WD distribution. 



tion p{E) = p\[E) for the same values 
of the coupling. 

To characterize the spectral statistics, we use the parameter rj. A priori, it is a function 
of the position E in the spectrum, of the coupling strength T and of the number m of 
intruders: ri{E / A,r / A,m). When F ^ A, the resonances do not overlap and is a 
single function of E/T as in the case of a single intruder (Fig. ^). When the resonances 
overlap enough (F > A), r] becomes almost energy independent: 77(F/A,m) (Fig. 
For infinite coupling, the statistics only depends on the number of intruders: 77(771). 
As the number m of these intruders increases, a smooth transition towards the WD 
regime is now observed. Fig.(^ shows how the parameter 77 depends on m and in the 
inset, it is shown that for m = 20, P{s) is indeed very close to the Wigner surmise. 
This smooth transition can also be observed in the function S]^(i?). (fig. |^). 




Fig. 5. E^(i5) for a Poissonian back- 
ground with an infinite coupling to 
m — 1,2,4,10,20 and 50 intruders, for a 
Gaussian coupling, F S> A. It approaches 
the WD curve when m increases. 
We have also calculated the distribution 
we compute the probability f{x) for the ^ 



Fig. 6. Function f{x) for a Poissonian 
background coupled to m = 1, 5 or 20 in- 
truders. As m increases, it approaches the 
Porter- Thomas distribution 

r the coefficients of the strength function: 
.ght Cn\ — |(A|n)p of an intruder |A) on an 
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eigenstate \n) to be equal to x. The random matrix theory predicts that f{x), when 
properly renormalized should be equal to f{x) = l/\/27rxe~^/^ (for a time reversal 
symmetric system) [|l9|. We have checked, that the function obtained is the same 
for an infinite coupling, or for a finite coupling when the eigenstates are taken in the 
central part of the Lorentzian. 

Consider now the situation where the number of intruders is large but the coupling is 
finite m ^ F/A. The statistics then depends only on the overlap between resonances: 
77(r/A) as it is shown on Fig.|^. Moreover the statistics of exact eigenstates do not 
depend on the relative position of the intruders: we have checked that if the intruders 
have all the same energy, or if they are randomly distributed, or if they are regularly 
and equally spaced, the statistics of eigenstates is always the same. 
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. A/A, = 5. 
• A/A, ^ 8. 
■ A/A, = 10. 
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Fig. 7. 77 vs r/A for different distances Fig. 8. P{s) for an infinite complex 

A between 10 intruders equally spaced. coupling, as the number of intruders in- 
The coupling is Gaussian. The horizon- creases, 
tal dashed line shows the value of rj in the 
limit of infinite coupling for 10 intruders. 

Finally, we have studied how the transition occurs when the coupling coefficients V\k 
are complex, which corresponds to a time reversal breaking situation. The modulus 
of V\k is chosen to follow a Gaussian distribution, and its phase follows a uniform 
distribution. As above we have considered the case of infinitely strong coupling. As in 
the time reversal symmetric case, the transition depends on the number of intruders. 
Eq. shows that the case of one intruder is identical to the time reversal symmetric 
case because the phases cancel. In particular the nearest spacing distribution is linear 
for small spacings. On the contrary, when there are two or more intruders, the spacing 
distribution is quadratic for small spacings, because the off-diagonal terms in eq. 
are now complex. As the number of intruders is increased, the distribution gradually 
evolves towards the distribution given by a random Gaussian unitary ensemble. 
To conclude, we stress that the cross-over to WD statistics is induced by indirect terms 
Mij, with i ^ j, in the matrix (Q). These off-diagonal terms of the form Vx-kVk\j are 
missing in a Cayley tree representation of the hierarchical structure considered in rcf. 
[0. They are however essential for a correct description of the spectral correlations. 



IV. Spectral function in the case of several intruders 

Up to now, we have only considered the statistics of the eigenvalues. We now consider 
the LDOS of the intruders. The LDOS for a given intruder Ai, is given by the 
imaginary part of GaiAi- The matrix elements Mij in eq.(^ are replaced by their 
imaginary part (neglecting a shift of order A3 due to the real part): 



7 



My = iV., 



ij 



k 



(8) 



which, as long as the density of the background levels is much larger that the densit- 
of intruders, can be approximated by: My = —iVij = —iTT{V\.kVkXj)iy3 The eq.([ 
reads: 



E - + ir. 



f^i^A fA 





V Glx,„ J 



whose solution is: 



{Xi\a){a\Xi) 

E — Er, 



(9) 



(10) 



The complex eigenvalues Ea give the position and the width of the resonances. We 
now consider several cases: 

1) When the coupling to intruders is uncorrelated and symmetric so that {VxikVkXj) = 
{Vxik){VkXj) = for i 7^ j, the off-diagonal elements in eq.(^ vanish and one finds 
that each level coupled to the continuum is broadened into a Lorentzian according to 
the Fermi golden rule; neither its center nor its width are altered: "the resonances do 
not talk to each other" . 

2) When the couphng is such that {Vx-kVkXj) = (Vxifc) (VfcAj ) 7^ 0, the off-diagonal 
elements of the matrix M are now finite and the resonances are coupled. Given the 
form of the off-diagonal matrix elements, the eigenvalues can be calculated easily and 
found to be solutions of: 



7j 



^ E^-Ex^+t{T,,-j,) 



(11) 



^jj ^ i't^3(^Ai*:V/cA ) are the width of the uncoupled resonances and jj = 
nMVx,k){VkxS- 

3) Consider first the very specific case where Tjj — jj . This case corresponds to taking 
a constant coupling Vx k — Vj (which may still depend on the intruder j). It has been 
studied recently by Konig et al. for the case of two intruders, e.g. two discrete levels 
of a quantum dot coupled to an electron reservoir pO| . The eigenvalues are solutions 



of eq. ( |11[ ) with Vjj = 7^ When the 7^ increase and become large compared to 
the distance between intruders A, the resonances are split into several peaks whose 
centers are given by: ^ Ex^ ~ and are placed at intermediate positions 

between the unperturbed levels. Thus the LDOS px^ consists in a series of (m — 1) 
peaks centered on the positions SR-Bq and whose widths are given by (^ / {Ex^ — 



Er, 



A~/7. The corresponding spectral weights are: 



1 



{E^-ex,YY. 



(12) 



3 (eA,--Ec)2 



The (m— 1) peaks actually carry (1 — 1/m) of the spectral weight. The rest is carried 
by a very broad Lorentzian of width 7^ . This generalizes the case m = 2 studied 
in ref. |2^] where , at large coupling, a peak of spectral weight 1/2 appears in the 
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middle of the two original discrete levels. More generally, fig.(Q.a) shows that the 
spectral weight tends to concentrate at positions placed between the original states. 




Fig. 9. LDOS of the 3'''' out of 4 intrud- 
ers whose unperturbed energies are figured 
by strait lines. All Fjj's and ^j's are equal. 
Tj = 5, while 7j = (Golden Rule case 



Fig. 10. 77 vs U/Ac when the Fock 
space is restricted to the the last 2 gen- 
erations (2n* — 1) and (2n* + 1). r? is also 
a function of the DOS S. 



4 or 5. 



4) In the more general case where jj ^ Tjj, the spectral function is still the su- 
perposition of (m — 1) peaks, centered at the same positions, but whose widths are 
Max{{Tjj — 7j), A^/7j). If the fluctuations of the coupling distributions of the cou- 
plings V\-k are sufficiently large so that Tjj — 7^ ^ A, the peaks will transform into 
one single Lorentzian, whose width is however smaller than the width Tjj predicted 
by the Golden Rule (eq.|). 

5) In the most general case, {VxikVkx^) ^ {V\.k){VkXj), the matrix M has to be 
diagonalized numerically. In view of the preceding arguments we still expect the 
width of the resonances to be smaller than the bare resonance width. 



V. Coupling between generations of higher order 

We now consider the general case of a coupling between two generations (2n— 1) and 
(2n + 1) in the hierarchy of states considered in ref. 0]. The spacing between levels 
in a given generation is pl] , p^ : A2n+i — {2n)\n\{n + l)!A^"+^/e^". However, due 
to the two-body nature of the interaction, a level of generation (2n — 1) is coupled 
to a small number of states of the next generation. The distance between levels of 
generation (2n + 1) connected by the interaction is of order: Ac ~ nA^ ~ 4nA^/e^. 
The direct coupling between states of a given generation can be neglected since the 
distance between such states is A^ ~ A^/e which is always much larger than A3. We 
now determine the condition under which the many-body states constituted of these 
two generations obey WD statistics. If this condition is obeyed the many-body states 
obtained by successive couplings from generations (1) to {2n + I) will also obey a 
WD statistics. As noticed in ref. 0, the DOS in a generation (2n + 1) is much larger 
than in the previous generation. Due to the finite size of the system, this is true up 
to n* = ■y/e/2A. Then further generations are coupled to generations with smaller 
density of states, which certainly do not affect the appearance of WD correlations. 
We are thus especially interested in the two generations (2n* — 1) and {2n* + 1) which 
have the same level spacing S ~ Ae""*" . Typically the distance between connected 
states is A^ ~ A/n*^ ~ A(A/e)3/2 > S. 
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A priori three parameters are relevant, the inter-level spacing 5, the spacing between 
states connected by the interaction Ac and the typical matrix element of the interac- 
tion U. It has been argued that the mixing of the states and the cross-over to a WD 
statistics depend only on the ratio U/ Ac [Q-^jl2|. Then, the cross-over is expected to 
occur when this ratio is of order unity. Since U A/g where g is the dimensionless 
conductance Ec/A, this gives A/g ~ A{A/ef/'^ so that e ~ g^/^ as found originally 
by Jacquod and Shepelyanskii and recovered by Mirlin and Fyodorov Pp2| . 



However, fig. 10 shows that the crossover is not uniquely driven by the ratio U / A^ 



but it also depends on the density of states 5. For given values of U and Ac, when 
Ac/ 5 increases, the transition is faster. This may appear surprising since the density 
of non-zero coupling elements decreases. However, at the same time, the distance 
between levels decreases, so that mixing neighboring levels becomes more efficient. 
A level of the generation (2n* — 1) is typically coupled to a level of generation (2n* -1-1) 
at a distance Ac S> 5. The indirect coupling between two neighboring levels of 
the same generation (distant of 5) typically necessitates high orders in perturbation 
U{U /Ac)^ where p is of order y/KJS . The mixing of these neighboring levels and 
the transition to WD thus involve many such high order processes [|ll| . We have not 
yet succeeded in finding the correct criterion which involves both energy scales Ac 
and 5. 
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